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Abstract. An intense radio flare on the dMe star AD Leo, observed with the Effelsberg radio telescope and spectrally resolved 
in a band of 480 MHz centred at 4.85 GHz is analysed. A lower limit of the brightness temperature of the totally right handed 
polarized emission is estimated as Tb ~ 5 x 10^" K (with values Tb ^ 3 x 10^^ K considered to be more probable), which 
requires a coherent radio emission process. In the interpretation we favour fundamental plasma radiation by mildly relativistic 
electrons trapped in a hot and dense coronal loop above electron cyclotron maser emission. This leads to densities and magnetic 



field strengths in the radio source of n < 



and B ~ 800 G. Quasi-periodic pulsations during the decay phase of 



the event suggest a loop radius of r ~ 7 x 10 cm. A filamentary corona is implied in which the dense radio source is embedded 
in hot thin plasma with temperature T > 2 x 10^ K and density next < 10~^n. Runaway acceleration by sub-Dreicer electric 
fields in a magnetic loop is found to supply a sufficient number of energetic electrons. 
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1. Introduction 

Several nearby flare stars occasionally show intense radio 
outbursts. The most prominent events were observed from 
AD Leonis, a very active single M dwarf star (Gliese 388, 
dM3.5e, ^ 3.5 x lO^" cm, d = 4.85 pc = 1.55 x 10^^ cm). 
Its radio flares have been detected mainly in the decimetric 
range at various observing wavelengths around 20 cm (Lang et 
al. 



1983 



1989 



Lang & Willson 



1986; White et al. 1986; Jackson et 



Gudelet al. |1989| ; Bastian et al. |1990| ; Abada-Simon 



al. 

et al. |1994a[ |1997| ). The estimates of the brightness temperature 
of these emissions range from Tb ^ 10^° K to Tb 10^'' K, 
which all imply a coherent emission process. The circular po- 
larization is often high, but can vary from event to event be- 
tween negligible values and « 100 %. The emission can be dif- 
ferent or even absent at a neighbouring frequency only ^ 10 % 
apart. On several occasions, dynamic spectra have shown a 
wealth of fine structures in the time-frequency plane (pulsa- 
tions, sudden reductions, spikes, fast-drift bursts), largely sim- 
ilar to solar dynamic radio spectra. All these properties — high 
polarization degree, narrow bandwidth, and fine structures — 
require a coherent emission process. 

AD Leo flare observations in the microwave range, per- 
formed mainly near 6 cm, are less numerous and have pro- 



duced a less coherent picture so far. Possibly the observed 
events belong to two or three groups. (1) There are long-lasting 
(t ^ 1 hour), relatively weakly polarized bursts of low to mod- 
erate flux den sity, F pcak = 1.4-2 mJy (Gary et al. |1987 ; 
Jackson et al. 1989; Rodono et al. 1990). If a source radius 



equal to the stellar radius is assumed, these flux densities cor- 
respond to brightness temperatures Tb ~ (0.2-2) x 10^ K, 
consistent with incoherent emission by mildly relativistic elec- 
trons. (2) A more impulsive, completely left-handed polarized 
event of w 5 min duration and flux density F^^ak ~ 30 mJy 
was observed by Gary et al. ( 1987 ). Of similar nature may 



be the impulsi ve ons et phase of the gradual burst observed by 
Rodono et al. (1990), which had w 2 min duration, a flux den- 
sity Tpoak ~ 30 mJy, and was seen only at 2 cm, not at 6 or 
20 cm (no information on the polarization was provided). If 
the source diameter is estimated from the rise time using the 
Alfven velocity ua ~ 7 x lO^Wc/^p cms^^ (where the ratio 
of electron cyclotron and plasma frequencies is assumed here 
to be tJc/'jJp ^ 1 for the purpose of an estimate), one obtains 
Tb - 7 X 10^" K for the burst observed by Gary et al. (|1987|) 



and Tb ^ 3 X 10^ K for the burst observed by Rodono et al. 



(199C). At least the former of these bursts requires a coherent 
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emission mechanism. (3) Lecacheux et al. ( 1993 ) observed one 
event of w 400 mJy peak flux and ~ 4 s duration, which was 
considered likely to originate from AD Leo. The burst was re- 
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solved by an acousto-optical spectrograph operating in the 4.5- 
5 GHz range at the Arecibo telescope. It was extremely narrow- 
banded, clearly discernible only between 4.78 and 4.91 GHz, 
and only moderately polarized (« 20 %). The brightness tem- 
perature was estimated as Tb « 3 x 10^''(i?,/rsourco)^ K, 
where it is clear from the short duration that the source ra- 
dius must have been much smaller than the stellar radius, 
''source ^ R*', again a coherent emission mechanism is re- 
quired. 

The stellar origin of these radio bursts has been confirmed, 
in the majority of cases, by simultaneous observing with up to 
three widely separated telescopes, by imaging with the VLA, 
by eliminating interference based on its characteristic appear- 
ances in dynamic spectra, or (in case of single dishes) by using 
two feeds simultaneously with one feed in an "off" position. 
Also the similarity of the fine structures with solar radio fine 
structures suggests a stellar origin. 

In the 2-6 cm range, a transition from coherent emission at 
long wavelengths to incoherent emission at short wavelengths 
appears to occur in the flare emission of AD Leo, similar to the 
situation at other flare stars (Bastian 1990 ) and the Sun, per- 
haps shifted to somewhat shorter wavelengths at flare stars. It 
is understood as being due to three effects. First, the character- 
istic frequencies of coherent emission (suJc/2tt; s = 1, 2 for 
the cyclotron maser and suJp/2tt; s = 1, 2 for plasma emis- 
sion) are limited to ^ 10 GHz by the available field strengths 
and plasma densities in stellar coronae. Second, a general in- 
crease of the ratio Wc/usp with decreasing height in the corona, 
which makes cyclotron emission processes more important and 
plasma emission less important or quenches the latter com- 
pletely. Third, increasing free-free absorption of plasma emis- 
sion and of low-harmonic cyclotron maser emission at increas- 
ing plasma dens ity (im plied b y high er observing frequency i/), 
see, e.g., Dulk ( 1985 ), Benz (1984). However, also similar to 
the solar case, coherent emissions observed in this range, albeit 
being rare and possibly becoming progressively weak at de- 
creasing wavelength, are likely to originate closest to, or from 
wi thin, t he primary flare energy release volume (e.g., Benz et 
al. |2001| ). 

We have therefore begun a search for bursts from AD Leo 
in the 6 cm range using the Effelsberg radio telescope and a 
bandwidth of w 500 MHz centred near 4.85 GHz, which pro- 
vides a sensitivity limit of ~ 50 mJy in dynamic spectra. A 
well developed, intense burst was detected, which warrants a 
detailed study, presented in this paper. The burst had a peak 
flux of « 300 mJy, duration of ^ 1 min, a broadband spectrum 
that varied slightly during the event, and was w 100 % right- 
handed circularly polarized. The estimated brightness temper- 
ature, Tb > 5 X 10^" K, implies a coherent emission process. 
The temporal modulations of the decay phase resemble quasi- 
periodic pulsations. 

Simultaneous soft X-ray observations using ROSAT were 
not successful, since the burst occurred during source occulta- 
tion by the Earth. 

Our analysis of the event focuses on the nonthermal radio 
emission mechanism. There are two mechanisms that can lead 
to the required high brightness temperature: electron cyclotron 
maser (ECM) emission and plasma emission (e.g., Dulk 1985; 



Kuijpers 1989). Since gyroresonance absorption of ECM emis- 
sion by ambient thermal electrons presents a severe problem for 
temperatures and densities supposed to be typical of M dwarf 
coronae, we consider plasma emission in detail. The high po- 
larization degree requires that the plasma emission arises at 
the fundamental. We assume that a population of mildly rel- 
ativistic electrons is trapped in a coronal loop, develops a loss- 
cone distribution, and excites electrostatic upper-hybrid modes 
(analogous to the canonical source model of broadband so- 
lar plasma emission). Calculations of the brightness tempera- 
ture of fundamental and harmonic plasma radiation for typi- 
cal M dwarf coronal parameters confirm the possibility of fun- 
damental emission. This mechanism implies smaller magnetic 
field strengths and higher plasma densities than ECM emission. 
Electron acceleration by a sub-Dreicer DC electric field at the 
chromospheric footpoints of a coronal loop is briefly consid- 
ered and found to be consistent with the requirements on the 
number and energy of accelerated electrons in the observed ra- 
dio flare. Finally, a consideration of the pulsations supports the 
assumption that the radio source was formed by particle trap- 
ping in a loop. 

The observations are presented in Sect. ^ the ECM and 
plasma emission processes are discussed and source parame- 
ters are estimated in Sect.||, the calculation of the brightness 
temperature for fundamental and harmonic plasma emission is 
given in Sect. ^ and the particle acceleration is considered in 
Sect. ^. The origin of the pulsations is discussed in Sect. ^ and 
the results are summarized in Sect.^. 

2. Radio observations of tlie AD Leo flare at 
Effeisberg 

The radio observations of AD Leo at Effelsberg lasted a total of 
about 27 hours on the days May 17-19, 1997. During the first 
two days the meteorological conditions were not favourable so 
that possible bursts could not be identified with certainty, but a 
well developed burst with a duration of ^ 1 min was recorded 
on May 19, 1997, 18.945 UT 

The observations were performed with the new 8192- 
channel autocorrelation system (AK) and in parallel with the 
pulsar receiving equipment. The centre frequency of the mea- 
surements was 4.85 GHz. The total effective bandwidth was 
500 MHz for the pulsar mode and 480 MHz for the spectro- 
scopic mode. Both observing systems recorded the two total 
power signals (LHC and RHC) of the four channel receiving 
system. 

2.1. The pulsar mode 

Left and right handed circular polarized radio power was de- 
tected from the two feeds of the receiver, one being centred on 
the source and the other one in an "off" position. The same 
centre frequency of 4.85 GHz was employed and the detec- 
tors operated at an effective bandwidth of 500 MHz. Using fast 
(10 MHz) V/f converters, the four frequency encoded intensity 
channels provided the input to the Effelsberg Pulsar Observing 
System (EPOS). The V/f frequencies of the four channels, typ- 
ically about 4 MHz, were continuously measured in 1.0 ms in- 
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Fig. 1. Time profile of the AD Leo flare observed at the 
Effelsberg 100 m radiotelescope on May 19, 1997, 18.945 UT 
at 4.85 GHz. Right and left handed circularly polarized flux 
on-source and off-source are plotted. 

tervals. Hence a resolution of2.5x 10^'' in detected power and 
1 .0 ms in time were the limits of the pulsar data aquisition as 
it was operated. For a system temperature of 30 K and a tele- 
scope efficiency of 46 %, about 6 counts per sample correspond 
to the achievable sensitivity of 32 mJy. After each collection of 
1024 data samples, these, together with a /is timestamp, syn- 
chronized to the station clock, were recorded on disk. 

Figure |l] shows the time profile of the event. Here the data 
were binned down to a time resolution of 1 s to obtain a better 
signal to noise ratio. The contribution of the ground radiation 
has been removed by a parabolic fit to the data. The r.m.s. noise 
after binning is about 2 mJy. Some low level structures, espe- 
cially before 18.94 UT, are visible in all four channels and are 
very likely caused by weather effects. Only the feed centred on 
the source observed the strong burst of radio emission, which 
was also highly polarized. While the typical flux density of the 
radio burst is ^ 200 mJy in the right handed polarization, the 
signal in the left handed polarization channel is below the 3 
sigma level of 6 mJy. The polarization degree of the detected 
signal is, therefore, larger than 97 %. 

2.2. The spectroscopic mode 

The new AK was used to record the left handed and right 
handed polarization with a large bandwidth. Four frequency 
bands were recorded, each 160 MHz wide. The centre fre- 
quencies of the individual bands were -180 MHz, -60 MHz, 
H-60 MHz, and H-180 MHz with respect to 4850 MHz. Each 
band consisted of 128 channels that were 1.25 MHz wide. The 
integration time was set to 128 ms. The reduction of the data 
was done by calculating {S — B)/B (where S is the signal and 
B is the background flux), and the calibration used a flux den- 
sity of 7.5 Jy for 3C286. Where the sub-bands overlapped, the 
16 overlapping data points (20 MHz) were skipped because dif- 
ferent baseline characteristics at the end of the bands prohibited 
their individual comparison. Differences of the flux density be- 
tween the edges of adjacent bands were corrected by applying 



linear baselines in such a way that the corrected value became 
equal to the mean of the adjacent edges. The four bands were 
then merged to the effective bandwidth of 480 MHz. In most 
cases the systematic uncertainty resulting from this procedure 
is below 30 mJy, but values of up to 100 mJy were occasionally 
noted. It is assumed that there is no other flare emission outside 
the flare event visible in Fig. |l]. A linear baseline along the time 
axis was subtracted to remove the background emission (i.e., 
the flux level before and after the considered event was set to 
zero). Finally, a method of unsharp masking invented by Sofue 
& Reich (1979) was used to remove "scanning effects" along 
the time axis. The rm.s. in quiet regions of the spectral data 
was found to be « 50 mJy. 

The resulting dynamic spectrum is shown in Fig. ^ for the 
right handed polarization. Due to internal problems of the new 
AK, we lost the signal during certain periods. These time inter- 
vals have been cleared in Fig. |[ No signals were detected in left 
handed polarization. The right handed data show a clear sig- 
nal between 18.945 and 18.965 universal time. Between 18.952 
and 18.955 UT the largest differences between adjacent bands 
were detected (mean ^ 80 mJy, peak 100 mJy) at frequencies 
near 4630 MHz. The data in this time interval have, therefore, 
the lowest accuracy. The other time slices show systematic er- 
rors of only 30 mJy in addition to the statistical error of 
w 50 mJy. 



2.3. Spectral and temporal details 

The bottom panel of Fig. ^displays the flux obtained in spec- 
tral mode integrated over the whole bandwidth. The compar- 
ison with the time profile obtained from the pulsar data set, 
binned down to the same time resolution of 125 ms, shows ex- 
cellent agreement after 18.956 UT, but also some deviations in 
the time interval 18.952-18.955 UT. These deviations reflect 
the systematic errors discussed above. The upper right panel 
shows the average spectra in each of the four time intervals in- 
dicated. Some spectral evolution during the lifetime of the burst 
is evident. 

The burst showed an irregular pulsating structure during its 
decay phase (18.956-18.961 UT), visible in the dynamic spec- 
trum (Fig.0 top left) as well as in the broad-band time pro- 
files (Fig.pTbottom). This structure is shown in greater detail 
m Fig.|.Inthisfi gure, we have subtracted the underlying flux 
level, obtained by taking the average over the whole band of the 
spectral data set at 18.956 UT and 18.9602 UT and linearly in- 
terpolating in between. An irregular sequence of about 6 pulses 
is clearly discernible in the subtracted time profile shown in the 
bottom panel. The subtracted spectrum (Fig. ^ top) reveals also 
an irregular spectral structure. The first three pulses do not ex- 
tend over the whole observed band and definitely show the flux 
maximum at different frequencies. The final two pulses seem 
to extend over the whole band, i.e., their bandwidth exceeds 
10%. 

The power spectrum of the pulse sequence (Fourier trans- 
form of the autocorrelation function of the background- 
substracted time profile of Fig.||) is shown in Fig.Q. It reveals 
a pulsation period of « 2 s (0.5 Hz), which is clearly different 



4 A.V. Stepanov et al.: Microwave plasma emission of a flare on AD Leo 



AD Leo 19-MAY-1997 (Effelsberg) 




Time Profile 

E 



18.945 18.950 18.955 18.960 18.965 
t [UT] 

Fig. 2. Plot of AK90 spectroscopic observations of the AD Leo flare on May 19, 1997 (top left). Intervals with no reliable spectra 
were cleared. The time profile averaged over the spectrum is shown in the bottom panel together with the pulsar (EPOS) data 
(binned down to the same time resolution and offset by 100 mJy for clarity). The top right panel shows average spectra in the 
four time intervals indicated. 




from the telescope oscillation period of « Is. Furthermore, 
since the pulses possess spectral structure, they cannot be of 
purely instrumental origin. 



3. Emission meclianism and source parameters 

3. 1. Characteristics of M dwarf coronae 

First we recall a few basic findings about the physical condi- 
tions in the atmospheres and coronae of M dwarf stars. Such 
stars possess outer convections zones and show universal X- 
ray emission (Schmitt 1997). Model fitting of low resolution 
X-ray spectra (Schmitt et al. |1990| ; Schmitt |1997| ) as well as 
studies of selected stars at significantly higher spectral res- 



olution (Hawley et al. 19951 has shown that the coronae of 



these stars contain plasma at a variety of temperatures with the 
emission measure usually rising to beyond T > 2 x 10^ K. 
Unfortunately, direct measurements of density in the coronae of 



M dwarf stars have not been made yet; Fossi & Landini (1994 ) 
quote an upp er limit of ^ 3 x 10^^ cm^'^ for AU Mic. Katsova 



et al. ( |1987| ) estimated the average density at the base of the 
corona from the X-ray emission measure of AD Leo, assuming 
a hydrostatic stratification and a temperature of T = 3 x 10^ K, 
as n ~ 5 X 10^ cm^^, a value probably smaller than the true 
value of density. From kinetic pressure balance, one expects the 
hot component to possess a significantly smaller average den- 
sity {n < 10^ cm^"^) and because of the high temperature the 
corresponding scale height is of order of a stellar radius. The 
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Fig. 3. Background-subtracted dynamic spectrum and time pro- 
file of the pulsating structure contained in the burst, obtained by 
integrating the subtracted dynamic spectrum over frequency. 

photospheric magnetic field strength of AD Leo is B ^ 3800 G 



with a high filling factor of order 75 % (Saar & Linsky 1985 ). 
The magnetic scale height is unknown; it is often assumed to 
be smaller than the density scale height, but such estimates are 
very model-dependent (see Saar & Linsky 1985 1. Similar to the 
Sun, the plasma-beta of the corona must be very small, even for 
B ~ 500 G we find = 4 x 10"'' < 1 for the density and 
temperature given above. Hence, already a weak twist in the 
magnetic field is able to support large perpendicular pressure 
gradients, for example, loops with T <; 2 x 10^ K and much 
higher densities of n ~ 5 x 10^^ cm"'^ embedded in an av- 
erage corona. A filamentary corona, where magnetically con- 
fined dense loops, generally cooler than the average due to en- 
hanced radiative or conductive losses, are embedded in a dilute 
hot corona of lower field strength is consistent with the obser- 
vation (White et al. 1994 ) that the hot component in dMe stars 
is permeated by relatively weak magnetic fields, B < 1000 G. 
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Fig. 4. Power spectrum of the time profile plotted in Fig.||. 

3.2. Estimate of the brightness temperature 

The brightness temperature of fully polarized radio emission is 
related to the flux density by 

"2 



(1) 



where n is Boltzmann's constant, d is the distance to the object, 
and A is the projected source area. Taking A ^ ttRI, one ob- 
tains Tb ^ 5 X 10^° K at peak flux. This is a very conservative 
lower limit on the brightness temperature, because the emis- 
sion is expected to arise from magnetically dominated struc- 
tures, e.g., a loop or a group of loops, which do not cover the 
whole stellar surface. 

By estimating the source extent from the timescale of the 
flux rise and a characteristic velocity, one can attempt to obtain 
a tighter limit on the brightness temperature. Let us first con- 
sider the velocity of the radiating mildly relativistic electrons, 
vi, which is appropriate if the expansion of a cloud of fast elec- 
trons in a more or less stable configuration, e.g. a magnetic 
trap, is seen as rise of the radio flux. The flare shows two main 
phases of flux rise. If extrapolated backwards to zero flux level, 
both phases last '-^ 7 s. Using vi ^ c/3 for ~ 50 keV particles, 
we obtain a source extent at peak flux of Lsourco ^ 7 x 10^° cm 
— about two stellar radii. As reviewed in Sect. [3.l| , one can- 
not expect that the high densities implied by the plasma emis- 
sion mechanism or the high field strengths implied by the ECM 
emission mechanism exist on such spatial scales; rather, loops 
of much smaller height are to be expected. Typically, loop 
heights are smaller than the pressure scale height (Rosner et al. 
19^, which is of order ~ i?, /3 for T - 3 x 10^ K. Therefore, 



the rise of the radio flux cannot be explained as a source expan- 
sion with the velocity vi. 

Alfvenic (or sub-Alfvenic) source expansion can be ex- 
pected in case of progression of the energy release process to 
neighbouring loops. Using ujc/cop ~ 0.5 (in accordance with 
the interpretation of the burst as plasma emission), we obtain 
source sizes of A^^^ = Lsourco ^ 2 x 10^ cm and brightness 
temperatures of Tb ^ 3 x 10^^ K. 

Alternatively, the rise of the flux is due to an increasing 
number of fast particles in the trap. For this case we suppose 
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a loop length £ ~ O.li?* based on the expectation that dense 
loops in general are cooler and possess a smaller pressure scale 
height than the average corona. The timescale of the flux pulsa- 
tions, evident in Figs. || and ^ provides an estimate of the loop 
radius. If the pulsations are caused by radial MHD-oscillations 
of a coronal loop (as suggested in Sect. their period is given 
by tpuisc ^ noop/t'A (e.g., Aschwanden 1987), which yields 
the loop radius as rjoop ~ 7 x 10^ cm. The resulting brightness 
temperatures range from Tb 4 x 10^^ K (for A ^ 2rioop/3, 
taking into account that the trapped particles develop a loss 
cone only in a section of the loop) to Tb ^ 10^^ K (for a source 
viewed along the loop axis, A ^ T^rf^op^ 

It should be noted that even higher Tb result in the lat- 
ter case if the source is filamentary at scales ^ Hoop (i-C-, 
if "loops" are composed of fine "threads"), as recent EUV 
observations of the solar corona suggest (Reale & Peres 
200C). Higher brightness temperatures would also be obtained 
for Alfvenic source expansion if the expansion is only one- 
dimensional (e.g., along an arcade of loops as in large solar 
flares) or under the assumption of ECM emission, which re- 
quires uJc/ojp ^ 1. Although there remains a substantial uncer- 
tainty in the value of the brightness temperature, all estimates 
imply that the radiation is coherent. 

Two coherent mechanisms are generally considered rele- 
vant in the context of stellar (and solar) radio emission: elec- 
tron cyclotron maser (ECM) emission and plasma emission. 
Both mechanisms are intrinsically narrow-banded and require 
an inhomogeneous source to explain broadband spectra. The 
inhomogeneity arises naturally in a magnetic trap, usually sup- 
posed to be a magnetic loop, where energetic particles de- 
velop a loss-cone distribution functi on (e. g., Dulk 1985 ). As 
discussed in detail by Bastian et al. ( 1990 1, both mechanisms 
do have the potential to explain a high brightness temperature, a 
high degree of circular polarization, and strong time variability. 
However, the absorption problem presents a severe difficulty 
for both mechanisms, particularly for the ECM mechanism. 



3.3. Electron cyclotron maser emission 

For various models of a loss-cone distribution function of the 
energetic particles and a wide range of background plasma tem- 
peratures [T ^ (0.2-1) X lO*" K] it was found that one can 
expect growth of the electromagnetic modes due to the ECM 
instability only at the first cyclotron harmonic, slightly above 
Wc (Sharma & Vlahos |l984l Melrose et al. |1984- The ex- 



traordinary (x) mode grows fastest in the range ujp/ujc ^ 0.4 
and the ordinary (o) mode grows fastest in the range 0.4 ^ 
Wp/ojc ^ 1- Maximum growth occurs at large propagation 
angles, 6'm ^ 70°, for both modes. Second and higher har- 
monics dominate the ECM instability for uup/ujc <; 1, but 
these waves cannot grow in reality since the electrostatic and 
whistler mode instabilities attain even higher growth rates in 
this range and exhaust the free energy contained in the fast 
partic le distribution (Sharma & Vlahos 1984; Winglee et al. 



[Winglee 1985 1). Furthermore, the ECM instability of mildly 
relativistic (^ 30 keV) electrons at the second and higher har- 
monics is too weak to survive collisional damping for back- 
ground temperatures T ^ 1.5x 10^ K (supposing a density 
ratio of accelerated electrons [ni] to background electrons [n] 
of ni/n < 10^3; see Fig. 3 in Sharma & Vlahos |1984| ). 

The modes that grow at the first harmonic are strongly 
damped by thermal electrons as their ray path crosses the sec- 
ond, third, and fourth cyclotron harmonic layers for decreasing 
magnetic field strength outside of the source. This gyroreso- 
nance absorption increases with increasing background plasma 
temperature and prevents propagation of the primary growing 
ECM emission through the low harmo nic layers in hot coronae. 
Detailed calculations (Stepanov et al. 1995, 1999 1 for realistic 
coronal temperatures [T ~ (0.3-3) x 10^ K] and densities 
{n ^ 10^-10^^ cm^'^) have shown that there is no escape win- 
dow for the X mode at the second harmonic (s — 2), only a 
very narrow escape window (9 ^ 10°) for the x mode at s = 3 
and the o mode at s = 2, and only a moderately narrow es- 
cape window (6 < 30°) for the o mode at s = 3; also the 
fourth harmonic contributes significantly to the damping. This 
is illustrated in Fig. || for parameters relevant to the burst under 
consideration (where we have assumed that the magnetic field 
gradient scale. Lb ~ 10^ cm, lies between the plausible values 
of the loop radius and the loop length discussed in Sect. 3.2). 



The optical depth scales roughly as To,x oc Lbti/lo. A narrow 
escape window of the o mode at s = 3 and nearly perpendic- 
ular propagation is often discussed, but does not exist for our 
parameters and probably also not for conditions of AD Leo in 
general (note that in the figure we have taken the minimum 
temperature of the hot coronal component and a rather small 
value of Lb)- 

ECM emission can reach the escape window at near- 
parallel propagation only through a substantial amount of an- 
gle scattering. Such scattering, which tends to isotrop ize th e 
primary radiation pattern, is possible (Stepanov et al. 1999), 
but its efficiency is difficult to estimate. Considering the nar- 
row angular interval of the ECM instability and the narrow 
escape window, a reduction of the brightness temperature by 
about one order of magnitude results from geometrical con- 
siderations alone. Non-perfect efficiency of angle scattering is 
expected to lead to a further significant reduction of the bright- 
ness temperature. 

Assuming that nearly all of the energy contained in the 
fast particles can be transferred to the maser emission, max- 
imum intrinsic brightness temperatures of ECM radiation at 
- 5 GHz of - 10^9 K are obtained (Melrose & Dulk 11984. 



1988). (It should be noted that a somewhat different result 



dominant growth of the x mode at the second harmonic for 
ojp/ujc ^ 0.6 - was obtained for cold background plasma 



The saturation levels of the ECM instability obtained from 
quasilinear theory (Aschwanden 1990) or from the condition 
of phase locking between the particles and waves (Wentzel & 
Aschwanden 1991) are strongly dependent on the density and 
energy of the energetic particles but reach values of ~ 10^^ K 
or higher for typical parameters. With these values, only a mod- 
erate efficiency of angle scattering (10^"^ . . . 10^^) is required 
for ECM emission at the estimated brightness temperature of 
^ 2 X 10^^ K (for Wc/^p ^ 2). Hence, one cannot exclude 
first-harmonic (s = 1) o-mode ECM emission, but the difficul- 
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Gyroresonance Absorption Mode Radiation Pattern 
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Fig. 5. Optical depths Ts for gyroabsorption of microwave 
emission of frequency i/ = uj/2tt = 4.85 GHz at differ- 
ent harmonics s vs. angle between the magnetic field and 
the wave vector of the ordinary and extraordinary mode for a 
Maxwellian plasma with T = 2 x 10'' K, n = 10^ cm^'^, and 
a magnetic field scale height Lb ^ cm. The uppermost 
lines refer to s = 2, and higher harmonics show monotonically 
decreasing optical depth. 



ties with the gyroresonance absorption do not suggest it to be 
the prime candidate as was often assumed in the literature. 

3.4. Plasma emission 

Due to the strong gyroresonance absorption of ECM, we con- 
sider only the plasma emission mechanism in the follow- 
ing, showing that it can, in fact, explain the observations in 
spite of the high observing frequency. For a trapped energetic 
particle distribution, plasma emission is based on the elec- 
trostatic upper-hybrid instability. This instability was investi- 
gated over a broad parameter range and for several slightly 
different loss-cone distribution functions of the fast particles 
by Zheleznyako v & Z lotnik ( |1975| ), Winglee & Dulk ( |1986| ), 
Stepanov et al. (1999), and many others (see the references 
therein). Only the range utp > ojc is relevant, since in the op- 
posite case, the growth rate of the electrostatic loss-cone in- 
stability drops significantly and the gyroresonance absorption 
comes into play. The electrostatic waves are always excited 
close to the upper-hybrid frequency based on the density of 
the background plasma, o^uh = i^p + ^^c)^^^' ^'^ wavenum- 
bers scaling with the cyclotron radius of the energetic parti- 
cles, kvi/ujc oJp/i^c, and at near-perpendicular propagation 
angles, 9^ ^ 80°. Nonlinear scattering processes lead to ra- 
diation at the fundamental, ujf « ujuh, and at the harmonic. 



"p/"c= 2 



/ cl/2 




Fig. 6. Angular pattern of fundamental o mode plasma emis- 
sion resulting from Rayleigh scattering of perpendicularly 
propagating plasma waves, according to Eq. (A.l). Radiated 
power is plotted at a linear scale. The arrow indicates the di- 
rection of the magnetic field. 

ti>h ~ 2(jJuh (see Sect. Q). We have oj^h ^ for all relevant 
values of the parameter LOp/uc- 

If thermal effects are taken into account, the electro- 
static waves and the fundamental emission occur at Wf « 
Wuh + 3fc^Wy/2(jJuh, where k is the wavenumber and vt = 
the electron thermal velocity. This frequency 
lies below the x mode cutoff frequency for LUp/ujc SS 20 (e.g., 
Bastian et al. 1990 1, which naturally leads to complete o mode 
polarization. Also harmonic plasma emission (which is much 
less affected by gyroresonance and free-free absorption than 
fundamental plasma emission) can lead to strongly polarized 
radio waves from hot coronae in the range 1 < (t^p/cjc)^ S5 5, 
but polarization degrees of k, 100 % over a considerable band- 
width, as observed in t he pre sent event, would be difficult to 
explain (Stepanov et al. 1999). Therefore, fundamental plasma 
emission has to be considered. 

Fundamental plasma emission at high frequencies suffers 
from free-free absorption due to the implied high plasma den- 
sities. The free-free optical thickness for fundamental plasma 
emission of frequency v — uOp/2-K (neglecti ng for the moment 
the influence of the magnetic field) is (DuUc 1985 ) 



Tff 



1.5 X IQ-^'^T-^/'^v'^Lr. 



(2) 



where the density in the source is given by u)p and Ln is the 
scale length of the density decrease along the ray path. For 
V — h GHz (n = 3 X 10^^ cm~^), the scale length leading to 
Tff = 1 becomes i„ = 2.7 x 10^"^ T"^/^ cm. It is immediately 
seen that the radiation could not escape if the radio source were 
surrounded by plasma at ~ 3 x 10^ K; the ray path must run 
primarily through the hot coronal component. It is also clear 
that even the hot component of M dwarf coronae (T < lO'* K) 
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would not permit escape of the radiation if the impHed high 
densities would slowly decrease along the ray path according 
to a hydrostatic atmosphere: Z„(r — 10® K) ^ 3 x 10^ cm. 
However, as noted above, such densities, being two orders of 
magnitude above the average, can exist in the corona of AD Leo 
only in thin filaments. If the external density is < 2 x 10^ cm^'^, 
we obtain Z„(T = 2 x 10^ K) > 3.5 x 10^° cm = per- 
mitting escape of the radiation. Since the average density of 
the hot coronal component is lower than this value and since 
refraction into regions of low density tends to "guide" the radi- 
ation through the hot component, we expect that fundamental 
plasma emission at 5 GHz can escape. 

Whether free-free absorption is important also within the 
source remains open. For i/ = 5 GHz, T — 2 x 10^ K, and 
Ln = Hoop ~ 7 X 10® cm, we have ts ~ 3. Under such 
conditions, the radiation can come from a surface layer of the 
source only. However, Tff < 1 for somewhat higher tempera- 
tures, T > 4 X 10^ K, which may exist even under quiescent 
conditions (Schmitt et al. 1990| ) or may be exceeded due to the 
flare heating. The calculations in the subsequent section sug- 
gest that the free-free absorption within the source looses im- 
portance at sufficiently high plasma wave energy densities (and 
brightness temperatures), where the emission is stimulated. 

Fundamental plasma emission can be influenced also by 
gyroresonance absorption. For Wc/wp ^ 0.5, the radiation has 
to cross the harmonic layer s — 3 and possibly even the layer 
s = 2. However, the effect is expected to be much weaker than 
for ECM emission in most cases, since fundamental plasma 
emission arises in the o mode, which is less damped than the 
X mode, and in two symmetric cones parallel and antiparallel 
to the magnetic field (Fig.^, i.e., in the direction of the es- 
cape windows. Depending on the actual value of oJc/uJp and on 
changes of the magnetic field direction outside of the inhomo- 
geneous source, the emission may escape without significant 
attenuation in some, perhaps many, cases and be completely 
absorbed in others. 

Finally, we remark that, also due to the high coronal tem- 
perature, the conversion efficiency of plasma waves into elec- 
tromagnetic radiation in a stellar environment exceeds the cor- 
responding solar effi ciency by one to two orders o f mag nitude 
(Abada-Simon et al. |19944 |l994bt Stepanov et al. |199^. 



3.5. Estimate of source parameters 

From the condition cJuIi/Stt — 4.85 GHz, and using again 
Wc/^p ^ 0.5 for the purpose of an estimate, we obtain the 
plasma density n ^ 2.3 x 10^^ cm~^ and the magnetic field 
B ~ 770 G. With T = 2 x 10^ K, the plasma-beta becomes 
/3 ~ 5 X 10~^ and the mean free path is Amfp ^ 10^ cm. With 
T = 5 X 10^ K, we have /3 - 0.13 and A,nfp 6 x 10^ cm. For 
both temperatures, the plasma is confined by the magnetic field 
and Amfp <C f source, hence the source can in fact be inhomoge- 
neous, as is required for broadband plasma emission. The range 
of densities and field strengths corresponding to the lower and 
upper limits of the observed band is 71 ~ (2.1-2.6)x 10^^ cm^-^ 
and B - 730-810 G. 



In comparison, if the burst were interpreted as ECM emis- 
sion, we would have uJc/^tt « 4.85 GHz, equivalent to 
B « 1700 G, and the condition ujp < ujc would lead to 
n < 2.9 X 10" cm-3 and (3 < 10"^ 

Minimum field strengths in the flare volumes of dMe stars 
have been derived from the observed X-ray luminosities and 
the estimated source volumes by van den Oord (1999). Values 
B ~ 500-2000 G appear to be typical. Both ECM and plasma 
emission are consistent with these estimates. 

4. Brightness temperature of plasma radiation 

In this section we consider the relative role of fundamen- 
tal (subscript f) and harmonic (subscript h) plasma radiation 
from a flaring loop on AD Leo. The plasma waves are ex- 
cited near the upper hybrid frequency, which is expected to be 
rather close to the plasma frequency since ojc < i^p- They are 
transformed into electromagnetic radiation at the fundamen- 
tal {ojf « LUuh ~ i^p) by the scattering on ions of the back- 
ground plasma (Rayleigh scattering). The conservation law for 
this scattering has the form 



LUf 



(fcf - k) ■ Vi 



(3) 



where uif, ui and fcf, k are the frequency and wave vector of 
electromagnetic and plasma waves, respectively, and Vi is the 
particle (ion) velocity. 

The second harmonic plasma radiation arises from a non- 
linear coupling of two plasma waves (Raman scattering), which 
obeys the resonance conditions 



k + ki = fch • 



(4) 
(5) 



The magnetic field is important in determining the fre- 
quency spectrum of the radiation which results from the con- 
version of plasma waves into electromagnetic waves, but it 
does not affect the total power of the radiation (Akhiezer et 
al. 1975). Moreover, for the loss-cone instability, the spectrum 
of plasma waves can be considered as isotropic (Zaitsev & 
Stepanov 1983; Stepanov et al. 1999). Hence, to estimate the 



brightness temperature, we can use the formulas for the weak 



magnetic field case given by Zheleznyakov (1996) (see also 



Zaitsev & Stepanov 1983) 



The solution of the transfer equation for the brightness tem- 
perature of radiation can be represented in the form 



Tb 



a 



1 — exp 



(A^c + Mnl) 



(6) 



Mc + Mm 

with the emissivities for the fundamental and the second har- 
monic 



36 v„ 



■w, 



tth 



(2^) 



5 ^3 ^2 



15^/3w>l 



nT 



(7) 
(8) 



and the absorption coefficients due to the electron-ion colli- 
sions 



fJ-c 



(9) 
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due to the scattering of the waves on background plasma ions B = 2 v 3 — (k 

c 



(10) c 



and due to the decay of an electromagnetic wave into two 
plasma waves 



Mnlh 



Here 



(27r)^ w 
15V3 vi £, 



5.5nr-3/2in(l04T3/2n-V3) 



(11) 



(12) 



is the electron-ion collision frequency, Wg = ^/2,k{l)vTc/ui-p 
and Wg — V3c/2 are the group velocities of the electromag- 
netic waves at the fundamental and at the second harmonic, 
respectively, vt — (kT/to)^/^ is the thermal velocity of the 
background electron component, vi = c[l — {m(? /{kTi + 
rnc^))^]^/^ is the velocity of the energetic electrons, m/M is 
the mass ratio, and w = W/nnT is the normalized energy 
density of the electrostatic waves, where W = J Wkdk. The 
formal parameter ^ is related to the width of the electrostatic 
wave spectrum by 



{Akf 



47r 



(13) 



We shall use fcmin ~ ^p/c and fcmax ~ ^p/5vt (which takes 
damping of the electrostatic waves at the thermal background 
into account and implies Ti ^ 25T). 

For dominance of fundamental emission, we require stim- 
ulated emission (maser effect), i.e., a negative absorption co- 
efficient. This happens for d{kWk) / dk > 0, which is valid 
if the plasma wave spectrum is sufficiently flat or has a pos- 
itive slope. This condition is easily satisfied for plasma wave 
turbulence spectra formed by nonlinear processes (Kaplan & 
Tsytovich [T973l ). 

The exponent in Eq. (^) is evaluated as r = J^^ /id/, as- 
suming a flat plasma wave spectrum. Since the density in the 
source is inhomogeneous with a scale length i„ — |n/Vn|, 
the integral is carried out only in a thin layer of depth 



3£. 



(14) 



in which the frequency of the electromagnetic waves is approx- 
imately constant. Since kf^^^ <C Loi/v^, we have L ^ L„. 
Furthermore, the relation dl = 6LnV^kdk / ujp is used (see Zai- 



tsev & Stepanov [ 1983 ] for details). The resulting brightness 
temperatures for fundamental and harmonic radiation are 



i^pi — Cw V 



Dwye 

l/ei + Fw/^ 



(15) 
(16) 



where we have introduced the abbreviations 



A 



TrUJn^ , 

36z;2 P ' 



D = 



F = 



TT m Vi 

2(27r)5 



nc 



15 

(2V3c)"^ 
2(27r)2 c 
15 



(18) 
(19) 

(20) 
(21) 
(22) 



The brightness temperatures are plotted in Fig. ^ for plasma 
parameters supposed to be representative of a flaring loop on 
AD Leo (T = 2 X 10^ K, Ti = 5 X 10^ K, L„ = 10^ cm), 
i.e., we suppose that the dense loop is heated up by the flare to 
about the temperature of the hot coronal component. The ther- 
mal level of plasma wave fluctuations for these parameters is 
Wth {Gn^nr'^)^^ w 8 x 10~^° (where td — vx/i^p is the 
Debye radius). An exponential increase of the brightness tem- 
perature of fundamental radiation (stimulated emission) occurs 
for w ^ 3 X 10~^. The two curves intersect at w* « 3 x 10^^, 
and Ri 2x 10^^ K. Fundamental emission dominates above 



this level of plasma turbulence. Our estimates in Sect. 3.2 sug- 
gested a brightness temperature Tb ^ 3 x lO^'^ K, which is 
roughly in accordance with the calculated minimum fundamen- 
tal brightness temperature, given the uncertainty of the estimate 
and the approximate nature of Eqs. ( [l^ and ([1^. 

Fundamental radiation reaches brightness temperatures ex- 
ceeding the estimated lower limit of ~ 5 x 10^" K also in the 
range 3 x 10~^ w < w*, in which a brighter continuum 
around 10 GHz due to harmonic plasma radiation would have 
remained unobserved. Such a situation cannot be excluded but 
does not appear very probable, since it is not supported by solar 
radio observations. 

The brightness temperature of fundamental radiation in- 
creases rapidly with increasing w, reaching Tbf > 10^*^ K 
within a short interval of w values. Strong turbulence effects 
will then limit the brightness temperature to Tbt << 10^^ K, 
at which level the energy density in the electromagnetic waves 
becomes comparable to the thermal (background) plasma en- 
ergy density. Viewed reversely, in the regime of dominating 
fundamental emission, the plasma wave energy density w stays 
within a rather small range. 

We have to ask whether the plasma wave energy density 
at saturation of the electrostatic upper-hybrid instability of the 
trapped particles can reach the level w ^ 3x lO^"^. A summary 
of the diverging answers to thi s ques tion, found in the literature, 
was given by Stepanov et al. ( 1999 ). We mention here only the 
estimates quoted by Dulk( 1985), It; ^ 10^^, which are consis- 
tent with our interpretation of the burst as fundamental plasma 
emission, and the estimate by Zaitsev et al. (1997), who con- 
sidered the case of saturation by moderate quasilinear diffusion 
(with almost isotropic energetic electrons): 



ni Ti liiCT vi 



n T a J'ci/C 



(23) 



Here a is the mirror ratio of the trap and C is the loop length. 
^^^•^ Taking the same parameters as in Fig.]?] with cr ~ 3 and 
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Fig. 7. Fundamental and harmonic brightness temperature vs. 
plasma turbulence level according to Eqs. (|l5]) and (jl^) for 
Wp/27r = 4.85 GHz, T = 2 x 10^ K, Ti = 5 x 10^ K (43 keV), 
and Ln = 10^ cm. 




Fig. 8. Brightness temperature (solid) and plasma wave level 
(dashed) at the crossing point, Tbf (it;*) = Tbh(w*) (cf. Fig. 
in T-Ln parameter space for Wp/27r — 4.85 GHz and Ti = 
25r. The contour Unes are a factor 10^/^ apart. 



£ ~ 3 X 10^ cm, the latter estimate yields w ~ 0.23 ni/n. 
Dominance of fundamental plasma emission then requires a 
density ratio ni/n ~ 10^''^, which is quite high but still rea- 
sonable (it can be reduced somewhat if Ti/T is larger). We 
conclude that the 4.85 GHz emission in the considered radio 
flare can be due to fundamental plasma radiation. 

A comparison with the case of decimetric emission is of in- 
terest. Repeating the calculation for a density that corresponds 
to a plasma frequency uj^j^-K — 1.4 GHz with the same values 
of T, Ti, and L„ gives the crossing point of the fundamen- 
tal and harmonic brightness temperatures at w* ^ 10^* and 

~ 3 X 10^^ K. Hence, harmonic plasma emission tends to 



dominate at longer wavelengths (Stepanov et al. 1999 ). 

An overview of the parameter dependence of the brightness 
temperature and plasma wave energy density at the intersec- 
tion Tbf(i«) = Tbh(w;) is given in Fig.^. It is seen that this 
brightness temperature increases with increasing background 
temperature and decreasing density gradient scale. Brightness 
temperatures Tb > 10^'^ K are achieved in a large range of 
parameter space. 

5. Electron acceleration 

In order to provide long-lasting radio radiation from a coronal 
loop, emitted by nonthermal electrons, a permanent source of 
energetic electrons is required, because the particles are lost 
from the magnetic trap due to diffusion into the loss cone. 
Under steady-state conditions, the acceleration rate N equals 
the loss rate into the loss cone. The particle loss rate due 
to pitch-angle diffusion of electrons by scattering on waves, 
which are excited by the loss-cone instability, is of the order 
for two-sided precipitation in a nearly symmetrical 
loop, where ripicc is the average number density of precipitat- 



ing electrons and S is the loop cross-sectional area at the foot- 
points. The fraction of precipitating electrons was estimated by 
Stepanov &Tsap (1999) as Tipicc ^ ^i/(2(t). Hence, 



N = niViS/a . 

Combining Eqs. ( ^3| ) and (^4|), we find 



In a ±1 



(24) 



(25) 



The level of plasma turbulence, w, can be estimated from the 
calculations of Sect. ^ Figure^ shows that w stays rather close 
to w* in the range where fundamental plasma emission domi- 
nates. 

We assume plasma parameters in the coronal loop as fol- 
lows: n = 2.3 X 10" cm-3, T = 2 x 10^ K, Ti = 5 x 10^ K 
(43 keV), £ = 3 X 10^ cm, S* - 10^^ cm^, cr = 3. From 
Fig. ^ we then have w* 3 x 10^^, and Eq. ( p5[ ) yields the ac- 
celeration rate required for the observed continuum brightness 
temperature as iV ^ 10^^ electrons s^^, which is nearly of the 
same order as the elect ron ac celeration rate during large solar 
flares (e.g.. Miller et al. |l997|). 



One possible origin for the particle acceleration in coro- 
nal magnetic loops on red dwarf stars is the production of run- 
away electrons in sub-Dreicer quasi-stationary electric fields, 
which can be a very efficient acceleration mechanism. Most 
authors who considered electron acceleration in DC electric 
fields have used the electric field of the currents in a plasma 
with the Spitz er resisti vity or with anomalo us resi stivity (e.g., 
Holman |l985t Tsuneta [l98^ . Zaitsev et al. ( |2000| ) considered 
another situation, namely the electron acceleration by electric 
fields that result from charge separation, which is caused (for 
strongly different electron and ion magnetization) by flows of 



A.V. Stepanov et al.: Microwave plasma emission of a flare on AD Leo 



11 



partly ionized plasma at the loop footpoints. They obtained the 
electric field component parallel to the magnetic field as 



Si 



1 — J- MfinVr Br 



^2 



B 



(26) 



where T = rinm^/ [11^171^ + nM) is the relative density of 
neutrals, is the "effective" frequency of ion-neutral colli- 
sions, Br is the small radial component of the magnetic field 
in a flux tube arising due to the expansion of the flux tube 
from the photosphere to the chromosphere, and Vr is the ra- 
dial component of the flow velocity. For sub-Dreicer fields, 
E\\ < Er. = e r^^ In A, the electron acceleration rate is 



N = 0.35n,M. 



Ed 



El 



3/8 



exp 



l2Eu 



El 



Ed 
4^ 



(27) 



(Knoepfel & Spong |1979| ). Here the index "a" denotes quanti- 
ties in the acceleration region, which has a volume Va- 

To estimate the resulting acceleration rate, we take the pa- 
rameters of the acceleration region in the stellar chromosphere 



as foUows: = 10^ K, 



10 



12 



and 



~ 5 X lO"* cms ^ (similar to fast solar convection). With 
Br 0.1 B we obtain Ei\ w lO^^ Vcm^i w Eu/'SO- An 
acceleration length h — 100 km is then sufficient for the ac- 
celeration of ~ 100 keV electrons. To obtain an acceleration 
rate of TV ~ 10^^ s~^, the cross section of the chromospheric 
acceleration volume has to be ~ 10^^ cm^, which is consis- 
tent with the chosen value of Br and the coronal loop radius 



estimated in Sect. 3.2 The density of energetic electrons re- 
quired for the observed radio fux was estimated in Sect. ^ as 
Til ~ 10~^n ~ 2 X 10^ cm^'^. This can be supplied by the 
chromospheric accelerator to a coronal loop with a volume of 
lO^^-lO^* cm^ in ~ 1 s. We conclude that this process is able to 
explain the acceleration underlying the observed intense radio 
burst. 

6. Origin of pulsations 

There are indications of flux pulsations with an average period 
of « 2 s during the decay phase of the event under considera- 
tion (Figs. ^ Usually four modulation mechanisms for stel- 
lar and solar radio emission are discussed in the literature: (1) 
modulation due to MHD oscillations of a flux tube containing 
the source, (2) intensity oscillations of the underlying electro- 
static waves during the nonlinear stage of the wave-wave or 
wave-particle interactions, (3) modulation of the particle accel- 
eration rate, and (4) modulation due to electric current oscilla- 



tions in an equivalent LRC-circuit (Aschwanden 1987 Zaitsev 
et al. 



1998) 



In order to obtain estimates of the radio source dimension. 



we considered the first mechanism of pulsations in Sect. 3.2, 
supposing radial oscillations of a loop which have a period 
ipuise ~ J'loop/wA, and obtained a quite reasonable value for 
the loop radius. Hoop ~ 7 x 10® cm. 

As for mechanism (2), the maximum value of the pulsation 
perio d can be esti mated as tn\ ~ 2Tr /lyei ~ 10^^ s (Zaitsev 



1971; Aschwanden 1987), which is too small in comparison to 



al 

3 X 



the observed period. 



The period of eigen-mode oscillations in a current-carrying 
coronal loop modeled as an equivalent LRC-circuit is iLC = 
(27r/c)(LC)^/2. The total inductance of a slender loop can be 
estimated as L w 10£ (Alfven & Carlquist 1967). The ca- 
pacitance depends on the current / flowing through the cross 
section of a loop, because just this current determines the self- 
consistent magnetic field and the effective dielectric permittiv- 
ity ofa loop plasma: C(/) w c'*nMS'2/(27r£/2) (Zaitsev et 
1998). Supposing again S ^ 10^^ cm^, we obtain <lc ~ 
10^^ / /. Taking iLC ~ 2 s, we find the value of the elec- 
tric current, I ^ 1.5 x lO^'^ A, and thus the magnetic field in 
the loop as B ^ 5 kG. This magnetic field value exceeds the 
photospheric field strength, which is highly improbable for the 
4.85 GHz source. 

Hence, only a modulation of the radio source by MHD 
oscillations or a modulated particle acceleration appear to be 
consistent with the observed period of the radio pulsations. It 
should be noted that MHD oscillations of a loop can lead also 
to a quasi-periodical regime of electron acceleration or electron 
injection into the loop. Alternative models exist, for example 
the quasi-periodic modulation of particle acceleration by dy- 
namic magnetic reconnection (Kliem et al. 2000). However, the 
fact that the pulsations are more clearly developed during the 
decay phase of the flare suggests MHD oscillations of a loop, 
which can be triggered by the primary flare energy release. 

7. Conclusions 

1 . Coherent plasma radiation can explain the intense broad- 
band AD Leo radio flare on 19 May 1997 at microwave fre- 
quencies (4.6-5.1 GHz), whose brightness temperature is esti- 
mated to be at least of order Tb ^ 5 x 10^° K (with values 
Tb <; 3 X lO^'^ K considered to be more probable). Electron 
cyclotron maser emission appears less probable than plasma 
emission, primarily due to the strong gyroresonance absorption 
in the hot stellar corona, but cannot be excluded. 

2. In case of the plasma emission mechanism, the ckcular 
polarization degreee of w 100 % requires emission at the fun- 
damental plasma frequency. Harmonic plasma radiation is less 
strongly polarized. 

3. The parameters of the radio source, assumed to be a 
magnetic loop, were estimated for plasma emission as n '--^ 
2.3 X 10" cm-^ and B - 800 G. Supposing T - 2 x lO'^ K, 
the plasma-beta is of order f3 ^ 5 x 10^^. 

4. A strongly inhomogeneous corona is required to avoid 
complete free-free absorption of the plasma emission. The ra- 
tio of source and external densities must be at least of order 
n/ucxt ~ 10^. 

5. Acceleration of electrons by DC electric fields at the 
footpoints of a loop, driven by photospheric convection, can 
give a rate of electron acceleration, N ~ 10'^^ s^^, sufficient 
to explain the observed flare in terms of plasma emission. 

6. A possible interpretation of the quasi-periodic modula- 
tion of the broadband radio emission during the decay phase of 
the flare is given by radial MHD oscillations of a flaring loop, 
which yields the loop radius as Hoop ~ 7 x 10® cm. 
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Appendix A: Angular pattern of fundamental 
plasma emission 

The angular pattern of fundamental plasma emission result- 
ing from Rayleigh scattering of perpendicularly propagating 
plasma wave s near the electron plasma frequency was obtained 
in Stepanov (197C) as 



p{e) = 



with 



C 



(A.l) 



A ^ ei-1- 

iV2 



B 



C 



1 



£2 



£l 



£3 



2 el 



sin^ e cos^ e 
{N^ sin^ 9 - £3)2 



iV2 



£1 



The elements of the cold-plasma dispersion tenso r, Sj, a nd the 
refractive index, N, are as given in Stepanov et al. ( 1999 ). Note 
that the assignment of the o and x modes to the sign of terms in 
the expression for N was erroneously exchanged in that paper 
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